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Abstract 

The World Summit on Sustainable Development (Johannesburg, 2002) encouraged the appli- 
cation of the ecosystem approach by 2010. However, at the same Summit, the signatory States 
undertook to restore and exploit their stocks at maximum sustainable yield (MSY), a concept 
and practice without ccosystemic dimension, since MSY is computed species by species, on the 
basis of a monospecific model. Acknowledging this gap, we propose a definition of "ecosystem 
viable yields" (EVY) as yields compatible i) with guaranteed biological safety levels for all time 
and ii) with an ecosystem dynamics. To the difference of MSY, this notion is not based on 
equilibrium, but on viability theory, which offers advantages for robustness. For a generic class 
of multispecies models with harvesting, we provide explicit expressions for the EVY. We apply 
our approach to the anchovy-hake couple in the Peruvian upwelling ecosystem. 

Key words: control theory; state constraints; viability; yields; ecosystem management; Peru- 
vian upwelling ecosystem. 

1 Introduction 

Following the World Summit on Sustainable Development (Johannesburg, 2002), the signatory 
States undertook to restore and exploit their stocks at maximum sustainable yield (MSY, see 
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(IClarkl. Il990f|). Th ough being criticized for decades, MSY remains a reference. Criticisms of MSY, 
like (jLarkinl . 119771 ). point out that MSY relies upon a single variable stock description (the species 
biomass), without age structure nor interactions with other species; what is more, computations 
are made at equilibrium. In fisheries, one of the more elaborate method of fixin g quotas, the ICES 
(International Council for the Exploration of the Sea) precautionary approach (jlCESl . 120041 ) , does 
not assume equilibrium (it projects abundances one year ahead) and assumes age structure; it 
remains however based on a monospecific dynamical model. Thus, in fisheries, yields are usually 
defined species by species. 



On the other hand, more and more emphasis is put on multispecies models (IHollowed. Bax. Beamish. Collie. Fog 



2000 ) and on ecosystem management. For instance, the World Summit on Sustainable Develop- 
ment encouraged the application of the ecosystem approach by 2010. Also , sustainabi l ity is a majo r 



goal of international agreements and guidelines to fisheries management (IFAO 



1999 



ICES 



20QJ). 



Our interest is in providing co nceptual insight as what could be sustainable yields for ecosys 



terns. In this, we follow the vein of (jKatz. Zabel. Harvey. Good, and Levin 



the c oncept of Ecologically Sustainable Yield (ESY), or of ( 



2003 ) which introduces 



Chapel. Deffuant. Martin, and Mullon 



20081 ) which defines yield policies in a 



viability approach. A general discussion on the 



approach to fisheries may be found in (jGarcia. Zerbi. Aliaume. Chi, and Lasserre . 



20031 ). 



ecosystem 



Our emphasis is on providing formal definition and practical methods to design and compute 
such yields. For this purpose, our approach is not based on equilibrium calculus, nor on intertem- 
poral discounted utility maximization but on the so-called viability theory, as follows. 

On the one hand, the ecosystem is described by a dynamical mo del controlled by harvest- 



ing. On the other hand, building upon (jBene. Doyen, and Gabay 



2001 



), constraints are imposed: 



catches are expected to be above given production minimal levels, and biomasses above safety bi- 
ological minimal levels. Sustainability is here defined as the property that such constraints can be 
maintained for all time by appropriate harvesting strategies. 



Such problems of dynamic control under c onstraints refer to viability (jAubinl . Il99ll ) or invariance 



(Oarke. Ledavev. Stern, and Wolenski 



19951 ) fra meworks, as well as to r 



or tubes for nonlinear discrete time dynamics in (jBertsekas and Rhodes! . 



eacha bility of target sets 



1971 



W e consider sustain a ble m a nagement issues formulated wit 



2001 



Bene and Doven 



lin such framework as in 



2003: 



Eisenack. Sheffran. and Kropp 



2006 
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Mullon. Curv. and Shannon 
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RaDaport, Terreaux, and Doven 
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De Lara. Doven. Guilbaud. and Rochet 
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A viable state is an initial condition for the ecosystem dynamical system such that appropriate 
harvesting rules may drive the system on a sustainable path by maintaining catches and biomasses 
above their respective production and biological minimal levels. We provide a way to characterize 
production minimal levels (yields) such that the present initial conditions are a viable state. These 
yields are sustainable in the sense that they can be indefinitely guaranteed, while making possible 
that the ecosystem remains in an ecologically viable zone; we coin them ecosystem viable yields 
(EVY). 

Thus, the EVY can be seen as an extension of the MSY concept in two directions: 1) from 
equilibrium to viability (more robust); 2) from monospecies to multispecies models. The second 
claim is obvious because, as we have recalled at the beginning, the MSY relies upon a single variable 
stock description (the species biomass), without interactions with other species. As for the first, 
recall that the MSY is the largest constant yield that can be taken from a single species stock over 
an indefinite period. By contrast, EVY are guaranteed yields, but they are not necessarily the 
annual catches. Indeed, it is by an adaptive catch policy (depending on the states of the stocks) 
that we shall be able to display yields indefinitely above the EVY. This is why we say that EVY 
are guaranteed yields, in the sense that catches cannot fall below the EVY. Viability can be seen 
as a robust extension of equilibrium: yields are not supposed to be sustained by applying fixed 
stationary catches, but are minimal levels which can be guaranteed by means of adaptive catch 
policies. 

The paper is organized as follows. In Section [21 we introduce generic harvested nonlinear ecosys- 
tem models, and we present how preservation and production constraints are modelled. Thanks 
to an explicit description of viable states, we are able to characterize sustainable yields. These 
latter are not defined species by species, but depend on the whole ecosystem dynamics and on all 
biological minimal levels. In Section [3l an illustration in ecosystem management and numerical 
applications are given for the hake-anchovy couple in the Peruvian upwelling ecosystem between 
the years 1971 and 1981. We conclude in Section [4] with possible extensions of the notion of ecosys- 
tem viable yields, on the one hand, to more general ecosystem models and, on the other hand, to 
bio-economic models so as to incorporate some economic considerations. We also discuss the limits 

3 



of the EVY concept. In the appendix, Section [A] is devoted to recalls on discrete-time viability and 
its possible use for sustainable management, while Section [B] contains the mathematical proofs. 



2 Ecosystem viable yields 

After a brief recall on the notion of maximum sustainable yield (for monospecific models), we 
introduce a class of generic harvested nonlinear ecosystem models, then present how to define 
maximum sustainable yields for this class. Next, we provide an explicit description of viable states, 
for which production and biological constraints can be guaranteed for all times under appropriate 
management. This makes possible to define ecosystem viable yields, compatible with biological 
and conservation constraints. We end up discussing relations between ecosystem viable yields and 
maximum sustainable yields. 

2.1 A brief recall on maximum sustainable yield 

We briefly sk etch the principles leadi ng to the notion of maxim um sustainable yield (MSY) (see 
199ol ) in continuous time and ( De Lara and Doyen . 



([Clark], 



200a) in discrete time). 



Consider a single population described by its total biomass B(t) at time t. Suppose that the 
time evolution of the biomass is given by a dynamical equation, a differential equation B(t) = 
Biol c (i?(t)) in continuous time or a difference eq uation B(t + 1) = Biol^(i?(t)) in discrete time. 



From this, build a Schaefer model (jSchaeferl . 



1954 ) by substracting a catch term h(t), giving B(t) 



Biol c (i?(t)) — h(t) or B{t + 1) = Biold{B(t) — h(t)). In general, to each biomass level B E (below 
the carrying capacity) corresponds a catch level h E = Sust(£? E ) for which the biomass B E is at 
equilibrium, solution of Biol c (£? E ) — h E = or Biol^(i? E — h E ) = B E . The maximum sustainable 
yield is the largest of such equilibrium catches: msy = max# E Sust(-B E ). 

2.2 Ecosystem biomass dynamical model 

For simplicity, we consider a model with two species, but it can be easiliy extended to N species in 
interaction. Each species is described by its biomass: the two-dimensional state vector (y, z) repre- 
sents the biomasses of both species. The two-dimensional control (v , w) comprises the harvesting 
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effort for each species, respectively. The catches are thus vy and wz (measured in biomass) [j] The 
discrete-time control system we consider is 




y(t)R y (y{t),z(t),v(t)) 
z{t)R z (y(t),z{t),w{t)) 



(1) 



where t stand for time (typically, periods are years), and where R y : R 3 — > R and R z : ]R 3 — > R are 
two functions representing growth factors (the growth rates being R y — 1 and R z — 1). This model 
is generic in that no explicit or analytic assumptions are made on how the growth factors R y and 
R z indeed depend upon both biomasses (y, z). 

In the above model, each species is harvested by a specific device: one species, one harvesting 
effort. This covers the multioutput settings case (e.g. several species in trophic interactions and 
targeted by the same fishing gear). Indeed, for this it suffices to state that both efforts are identical: 
v(t) = w(t) for all t = to, to + 1, 

2.3 Preservation and production sustainability 

We now propose to define sustainability as the ability to respect pr eservation and production 



minim al levels for all times, building upon the original approach of (jBene. Doven. and Gabav 
200ll ). Let us be given 



• on the one hand, minimal biomass levels B y > 0, B z > 0, one for each species, 

• on the other hand, minimal catch levels C y > 0, C z > 0, one for each species. 

A couple (jjo,zo) of initial biomasses is said to be a viable state if there exist appropriate 
harvesting efforts (controls) (v(t),w(t)) , t = to, to + 1, . . . such that the state path (y(t), z(t)) , t = 
to, to + 1, . . . starting from [y(to), z(to)) = (yo, zq) satisfies the following goals: 

• preservation (minimal biomass levels) 

biomasses: y(t) > B° y , z(t) > B\ , Vi = to, t + 1, . . . (2) 



1 In fact, any expression of the form c(y, v), instead of vy, would fit for the catches in the following Proposition [2] 
as soon as v \— > c(y,v) is strictly increasing and goes from to +oo when v goes from to +oo. The same holds for 
d{z, w) instead of wz. 
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and production requirements (minimal catch levels) 

catches: v(t)y{t) > & v , w(t)z(t) > & z , V* = to, t + 1, . . . (3) 



The set of all viable states is called the viability kernel (jAubinl . Il99ll ) . Characterizing viable states 
makes it possible to test whether or not minimal biomasses and catches can be guaranteed for all 
time. 

Here, by guaranteed, we mean that the yields have to indefinitely above the EVY, as reflected in 
the inequalities ([2]) and ([3]). We insist on the fact that, in the above definition of viable states, we 
say "there exist appropriate harvesting efforts (controls)": that is, to guarantee the EVY, we need 
to resort to adaptive catch policies (depending on the states of the stocks as in Corollary [6] in the 
Appendix). Viability can be seen as a robust extension of equilibrium: yields are not supposed to 
be sustained by applying fixed stationary catches, but are minimal levels which can be guaranteed 
by means of adaptive catch policies. 

Notice that, in the multioutput settings case, we need to add the constraint v(t) = w(t) for all 
t = to, t(y+l, . . .. Therefore, with an additional constraint, the set of viable states in the multioutput 
settings case will be smaller than the one considered above. 

The following definition summarizes useful and natural properties required for the growth factors 
in the ecosystem model. 

Definition 1 We say that growth factors R y and R^in the ecosystem model (TTJ) are nice if 
the function R y : R 3 — > M is continuously decreasing^ in the harvesting effort v and satisfies 
lim !) _ 5>+00 R y (y, z, v) < 0, and if R z : M 3 — > M is continuously decreasing in the harvesting effort w, 
and satisfies lim^^+oo R z (y, z, w) < 0. 

The following Proposition [2] gives an explicit description of the viable states, under some con- 
ditions on the minimal levels. Its proof is given in § IB.ll in the Appendix. 

The Proposition [2] may easily be extended to species in interaction as long as each species 
is harvested by a specific device: one species, one harvesting effort. However, it is not valid in the 
multioutput settings case. Indeed, it is crucial to have two distinct controls v(t) and w(t) for the 



2 In all that follows, a mapping (p : R — > M is said to be increasing if x > x' => ip(x) > ip(x'). The reverse holds for 
decreasing. Thus, with this definition, a constant mapping is both increasing and decreasing. 
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proof. Assuming that v(t) = w(t) for all t = to, to + 1, . . . would require other types of calculations 
for the viability kernel. This is out of the scope of this paper. 

Proposition 2 Assume that the growth factors in the ecosystem model ([1]) are nice. If the biomass 
minimal levels By, B\, and the catch minimal levels C y , C\ are such that the following growth factors 
values are greater than one 

R y {B y ,Bl% > 1 and R z (B y ,B z ,^) > 1 , (4) 

then the viable states are all the couples (y, z) of biomasses such that 

V>B y , z>B b z , yR y (y,z,%>B b y , zR z (y,z,^)>B b z . (5) 



Let us comment the assumptions of Proposition [2J That the growth factors are decreasing 
with respect to the harvesting effort is a natural assumption. Conditions (j4]) mean that, at the 
point (B y ,B z ) and applying efforts u = -gf-, v = -gf-, the growth factors are greater than one, 
hence both populations grow; hence, it could be thought that computing viable states is useless 
since everything looks fine. However, if all is fine at the point (B b ,B b z ), it is not obvious that this 
also goes for a larger domain. Indeed, the ecosystem dynamics given by (pQ) has no monotonocity 
properties that would allow to extend a result valid for a point to a whole domain. What is more, 
if continuous-time viability results mostly relies upon assumptions at the frontier of the constraints 
set, this is no longer true for discrete-time viability. 

2.4 Ecosystem viable yields 

Considering that minimal biomass conservation levels are given first (for prominent biological is- 
sues) , we shall now examine conditions for the existence of minimal catch levels 
First, we define (when they exist) the ecosystem viable yields. 

Definition 3 Let biomass conservation minimal levels B y > 0, B z > be given. Suppose that the 
growth factors in the ecosystem model ([T]) are nice, and that they take values greater than one in 
the absence of harvesting, namely: 

Ry(B*,B b z ,0) > 1 and R z {B y ,B z ,0) > 1 . (6) 
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Define equilibrium catches as the largest nonnegativ^ catches Cy*,C z '* such that 

Ry(BlX,^£) = 1 and R z (B b y ,B b z ,^) = 1 . (7) 

V * 

For a couple (yo,zo) of biomasses, define (when they exist) the ecosystem viable yields (EVY) 
Cy*{yo,z ) andC z '*(y ,z ) by 

Cy*(y , z Q ) := max{C y e [0, C y '*] \ y Ry{y , z , -yf) > B y } , 

(8) 

[ d'*(y , zo) := max{C z € [0, &/] \ z R z (y , z , > B b } . 
The term ecosystem viable yields is justified by the following Proposition UJ 

Proposition 4 Assume that the growth factors in the ecosystem model (pQ) are nice. For a couple 
(yo> z o) °f biomasses above preservation minimal levels - that is, yo > B y and zq > B\ - and 
satisfying 

y R y {yo,z ,0) > B b y and z R z (y , z , 0) > B\ , (9) 
the ecosystem viable yields C y '*(yo, zq) and C z '*(yo, zq) in (jHJ) are well defined. 



u 

< C y < Cy *(yo, Zo) and < C\ < C z '*(yo,Zo). Then, starting from the initial biomasses 
(y(to), z(to)) = (yo,zo), there exists appropriate harvesting paths which provide, for all time, at 
least the sustainable yields C y and C\ and which guarantee that biomass conservation minimal 
levels B y >0,B b z >0 are respected for all time. 

From the practical point of view, the upper quantities C y '*(yo, z$) and C z '*(yo, zq) in (jSJ) can- 
not be seen as catches targets, but rather as crisis limits. Indeed, the closer to them, the more 
vulnerable, since the initial point is close to the viability kernel boundary. 

Notice that the yield C y '*(yo, zq) depends, first, on both species biomasses (yo,Zo), second, on 
both conservation minimal levels B y and B\, third, on the ecosystem model by the growth factor 
R y ; the same holds for C z ' *(yo, zq). Thus, these yields are designed jointly on the basis of the whole 
ecosystem model and of all the conservation minimal levels; this is why we coined them ecosystem 
viable yields. 



What is more, consider catches C y and C z lower than these ecosystem viable yields, that is, 



3 Such catches are nonnegative because the growth factors in the ecosystem model (QJ are nice, hence continuously 
decreasing in the harvesting effort, and by <j6j) . 



s 



This observation may have practical consequences. Indeed, the catches guaranteed for one 
species depend not only on the biological minimal level of the same species, but on the other species. 
For instance, in the Peruvian upwelling ecosystem, it is customary to increase the biological minimal 
level of the anchovy before an El Nino event, but without explicitely considering to lower catches 
of other species. Our analysis stresses the point that minimal levels have to be designed globally 
to guarantee sustainability for the whole ecosystem. 

2.5 Ecosystem viable yields and maximum sustainable yields 

Now, we show how ecosystem viable yields are related to maximum sustainable yields. 

An equilibrium of the ecosystem model ([I]) is a couple (y E , z E ) of biomasses (state) and a couple 
(v E ,w E ) of harvesting efforts (control) satisfying 

y E = y E R y (y E ,z E ,v E ) , 
z E = z E R z (y E ,z E ,w E ) . 
The maximum sustainable yields, MSY y for species y and msy z for species z, are given by 

MSY„ := m&xv E y E and MSY^ := maxw E z; E . (11) 

v E ,w E v e ,w e 

They must be jointly defined because the ecosystem equilibrium equations (jlOD couple all variables. 

Say that the maximum sustainable yields MSYy and MSY^ are viable maximum sustainable yields 
if the corresponding biomasses equilibrium values y E and z E are such that y E > B y and z E > B\. 
In this case, MSYy and MSY 2 are ecosystem viable yields for the couple (y E , z E ) of initial biomasses: 
indeed, the stationary harvest strategy v(t) = v E and w(t) = w E drives the ecosystem model ([T]) at 
equilibrium (y E ,z E ) which satisfies the conservation minimal levels y E > B y and z E > B\. 

Notice that the maximum sustainable yields MSYy and MSY 2 are defined independently of the 
initial biomasses, whereas the ecosystem viable yields (EVY) Cy *(yo, zq) and ct'*(yo, zq) explicitely 
depend upon them. 

3 Numerical application to the hake— anchovy couple in the Peru- 
vian upwelling ecosystem (1971—1981) 

We provide a viability analysis of the hake-anchovy Peruvian fisheries between the years 1971 and 
1981. For this, we shall consider a discrete-time Lotka-Volterra model for the couple anchovy (prey 
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y) and hake (predator z), then provide an explicit description of viable states. 

We warn the reader that our emphasis is not on developing a "knowledge" biological model 
to "faithfully" describe the complexity of the Peruvian upwelling ecosystem. This formidable task 
is out of our competencies, and is not necessary for our analysis. Indeed, our approach makes 
use, not of "knowledge" models, but of "action" models; these are small, compact, models which 
capture essentials features of the system in what concern decision-making. In our case, we needed 
a compact model able to put in consistency biomass and catches yearly data between the years 



1971 and 1981. We chose a discrete-time Lotka-Vo 



candidate for a "knowledge" biological model (jHall 



terra model, despite well-known criticisms as 



1988 



Murray 



2002 ). but for its compactness 



qualities and for the reasonable fit (see Figure CQ) . 



3.1 Viable states and ecosystem sustainable yields for a Lotka— Volterra system 

Consider the following discrete-time Lotka- Volterra system of equations with density-dependence 
in the prey 

f y(t + l) 



R 

Ry{t) - -V 2 {t) - ay(t)z(t) - v(t)y(t) 



(12) 



k z(t + l) = Lz(t) + (3y{t)z(t) -w{t)z(t) , 

where R> 1, < L < 1, a > 0, /3 > and k = -^jK, with K > the carrying capacity for prey. 
In the dynamics ([I]), we identify R y (y, z,v) = R — —y — az — v and R z (y, w) = L + j3y — w. 
By Proposition U we obtain that, for any initial point (yo,Zo) such that 



yo>B b y , zo>Bl, y (R--y -az )>Bl, 



(13) 



the ecosystem sustainable yields are given by 

' Cl'*(y , z ) = min {B b y (R - £ B\ - olB\) - B y ,y (R - f y - az ) - B y ) 



cl'*(y ,zo) = Bl(L + PB V - 1) . 



(14) 



In other words, if viably managed, the ecosystem could produce at least Cy*(yo, zq) and C z '*(yo, zq) 
while respecting biological minimal levels B y and B\. 
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3.2 A viability analysis of the hake— anchovy Peruvian fisheries between the 
years 1971 and 1981 

The Peruvian upwelling ecosystem is extremely productive and dominated by anchovy (Engraulis 
ringens) d ynamics. It is well known that anchovy fisheries are very sensitive to envir onmental 



variability (jChecklev. Alheit. Qoseki. and Rovl . 



2009 



Sun. Chiang. Liu, and Chang 



2001 



), and the 



Peruvian anchovy is 



subject to environmental perturbations su c h as 



(ENSO) variability (jBertrand. Segura. Gutierrez, and Vasquea . 



11 Nino Southern Oscillation 



2004J) ) . However, for this simple 



predator-prey model using hake and anchovy, we have assumed that no uncertainties affect the 
ecosystem dynamics. Indeed, we feel that we have to go step by step in introducing the EVY 
concept, first focusing on the deterministic case. Thus, the period between the years 1971 and 
1981 is suitable for this first version of the model, due to the absence of strong El Nino events 
in the middle of the period. Furthermore, the long-term dynamics of the Peruvian upwelling 
ecosystem is dominat ed by shifts between altern ating anchovy and sardine regimes that restructure 



the entire ecosystem ([ Alheit and Niquen . 



20041 ). The period from 1970 to 1985 was characterized 



by p ositive temperature anom alies and low anchovy abundances, after the anchovy collapse in 



1971 ([Alheit and Niquen 



2004 ). so the competition between the fishery and hake was reduced due 



to low anchovy catches and anchovy mortality due to hake predation increase. Particularly, the 
changes in the ecosystem after 1971 led to an increase of five times, in average, in the predatio n 



rates of hake over anchovy between 1971 to 1980, with a peak in 1977 ([Pauly and Palomarea . 



19M)- 



Between the years 1971 and 1981, we have 11 couples of biomasses, and the same for catches. 
The 5 parameters of the Lotka-Volterra model are estimated minimizing a weighted residual squares 
sum function using a conjugate gradient method, with central derivatives. Estimated parameters 
and comparisons of observed and simulated biomasses are shown in Figure [TJ 

We consider values of B 9 = 7, 000, 000 t and B\ = 200, 000 t for minimal biomass levels 



(IMARPE 



2000. 



20041 ) . Conditions (113p are satisfied and the expressions (114p give the ecosystem 
viable yields (EVY) 

Cj'*(y , zo) = 5, 399, 000 t and &/(y , z ) = 56, 800 t . (15) 

In other words, such yields were theoreticaly susceptible to be guaranteed in a sustainable way 
starting from year 1971. In reality, the catches of year 1971 were very high and the biomasses 
trajectories were well below the biological minimal levels for fourteen years. 
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Observed 

Lot,ka-Volt,erra with 
density-dependence in prey 




- Observed 
• Lotka-Volterra with 
density-dependence in prey 




1970 1972 1974 1976 1978 1980 1982 



(a) Anchovy 



(b) Hake 



Figure 1: Comparison of observed and simulated biomasses of anchovy and hake using a Lotka- 
Volterra model with density-dependence in the prey (1971-1981). Model parameters are R = 



2.25 year" 1 , L = 0.945 year" 



67, 113 x 10 3 t (K = 37, 285 x 10 3 t), a = 1.220 x 10" 6 t" 



P = 4.845 



10" 8 t- 1 . 



The 4, 250, 000 t anchovy quo ta and the 55, 000 t hake quota, respectively, established for 



the year 2006 (jPBODUCE 



20051 ). or the 5,000 ,000 t anchovy q uota and the 35,000 t hake 



quota, respectively, established for the year 2007 (|PRODUCEl . 



20061 ) are rather close to the EVY 



cy*(yo 



figures C 



5,399,000 t and Cl'*(y ,z ) = 56,800 t. Thus, our approach provides reasonable 



4 Conclusion 



We have defined the notion of sustainable yields for ecosystem, and provided ways to compute 
them by means of a viability analysis of generic ecosystem models with harvesting. Our analysis 
stresses the point that yields should certainly be designed globally, and not species by species as 
in the current practice, to guarantee sustainability for the whole ecosystem. Our results have then 
been applied to a Lotka-Volterra model using the anchovy-hake couple in the Peruvian upwelling 



4 At this stage, we do not claim that the figures may be proposed as yields for the present management 
of hake-anchovy Peruvian fisheries. Indeed, our computations of EVY rely upon a dynamical model adjusted 
for some thirty years ago. To propose EVY, we should first dispose of a dynamical model adapted to the 
current situation, because it ought to reflect the new ecosyste m functioning and the depleted state of stocks 



( Ballon. Wo snitza-Mcndo. Gucvara-Carrasc o. and Bertrand , 



20081 ) . This is beyond the scope of this paper. 
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ecosystem. Despite simplicitjQ of the models considered, our approach has provided reasonable 
figures and new insights: it may be a mean of designing sustainable yields from an ecosystem point 
of view. 

We now discuss the limits of the EVY concept, as presented in this paper: application to biomass 
ecosystem models without age or spatial structure, no economic consideration, no uncertainties. 
We stress that EVY is a flexible concept, and we hint at possible extensions to incorporate the 
missing dimensions listed above. 

The framework we propose is not restricted to two populations, each described by its global 
biomass, but it may be adapted to several species, each described by a vector of abundances at 
age, or by vectors of abundances at age for each patch in a spatial model, etc. Suppose that the 
time evolution is given by a dynamical equation reflecting ecosystemic interactions and driven by 
efforts or by catches. Suppose that minimal safety levels (reference points) are fixed for biological 
indicators like spawning stock biomass, abundances at specific ages, etc. (such reference points for 
biological indicators like spawning stock biomass are generally given by international bodies like 
the ICES, or nationally). Ecosystem viable yields are minimal harvests for each species which can 
be guaranteed for all times while respecting the above minimal safety levels for biological indicators 
for all times too. 

It is often objected with reason that the MSY concept is developed without any economic 
consideration. As presented here, the EVY suffers the same criticism. However, the EVY concept is 
flexible enough to incorporate some economic considerations. For instance, upper bounds for fishing 
costs may be incorporated as constraints to be satisfied for all time, aside with minimal biomass 
levels. In this sense, EVY will be guaranteed yields compatible with biological and economic 
restrictions. 

As presented here, the EVY framework supposes that no uncertainties affect the ecosystem 
dy namics. Though we have the tools to tackle such an important issue (see stochasti c viability 



in ( De Lara and Doven 



2008; 



De Lara and Martinet! . 



2009 



Doven and De Lara 



20ld )). we feel 



that we have to go step by step. This paper introduces the EVY concept in the deterministic case, 



In addition to hake, there are other important predators of anchovy in the Peruvian upwelling ecosystem, 
such as mackerel and horse mackerel, seabirds and pinnipeds, which were not considered. Also, anchovy has 
been an important prey of hake, but other prey species have been found in the opportunistic diet of hake 
( Tarn. Purca. Duarte. Blaskovic. and Espinozal . 



2006) 
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providing an extension of the MSY concept in two directions: from equilibrium to viability (more 
robust), from monospecies to multispecies models. The extension to the uncertain case is currently 
under investigation. 

Thus, control and viability theory methods have allowed us to introduce ecosystem consider- 
ations, such as multispecies and multiobjectives, and have contributed to integrate the long term 
dynamics, which is generally not considered in conventional fishery management. 
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A Discrete— time viability 

Let us consider a nonlinear control system described in discrete-time by the dynamic equation 



where the state variable x{t) belongs to the finite dimensional state space X = M nx , the control 
variable u(t) is an element of the control set U = M. m while the dynamics f maps XxD into X. 

A controller or a decision maker describes "acceptable configurations of the system" through a 
set DcXxU termed the acceptable set 



' x(t + 1) = f(x(t),u(t)) foralliGN 
x(0) = xq given, 



(16) 



(x(t),u(t)) € D for all t G N 



(17) 



where B includes both system states and controls constraints. 
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The state constraints set Y° associated with B is obtained by projecting the acceptable set D 
onto the state space X: 

V° := Proj x (B) = {x G X | 3u G U , (x, u) G B} . (18) 

Viability is defined as the ability to choose, at each time step t G N, a control u(t) G U such that 
the system configuration remains acceptable. More precisely, viability occurs when the following 
set of initial states is not empty: 

' 3(m(0),ii(1),...) and (x(0), x(l), . . .) 



V(/,B) := 



Xq G 



satisfying 



The set V(/, B) is called the viability kernel (jAubin . 



(19) 

D and ([IT]) 

associated with the dynamics / and the 



1991 



acceptable set B. By definition, we have V(/,D) C V° = Proj x (B) but, in general, the inclusion is 
strict. For a decision maker or control designer, knowing the viability kernel has practical interest 
since it describes the initial states for which controls can be found that maintain the system in an 
acceptable configuration forever. However, computing this kernel is not an easy task in general. 

We now focus on some tools to achieve viability. A subset V is said to be weakly invariant for 
the dynamics / in the acceptable set B, or a viability domain of / in B, if 



VxGV, 3uGU, (x,u) G B and f(x,u) G V 



(20) 



That is, if one star ts from V, an acceptable control may transfer the state in V. Moreover, according 



to viability theory (jAubin 



1991 



), the viability kernel V(/, B) turns out to be the union of all viability 
domains, or also the largest viability domain: 

V(/, B) = U|V, V C V°, V viability domain for / in B | . (21) 

Viable controls are those controls u G U such that (x,u) G B and f(x,u) G V(/, B). 

A major interest of such a property lies in the fact that any viability domain for the dynamics 
/ in the acceptable set B provides a lower approximation of the viability kernel. An upper approx- 
imation Vfc of the viability kernel is given by the so called viability kernel until time k associated 
with f inB: 

' 3 (u(0),u(l),... ,u(fc)) and (x(0),x(l), . . . ,x(k)) 1 

satisfying (flgj) for t = 0. . . . . k - 1 \ ■ (22) 

and (H7J for t = 0, . . . , k 



¥ 



k ■- 



Xq G 
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We have 

V(/, B) C V fe+ i C Y k C V = V° for all k £ N . (23) 
It may be seen by induction that the decreasing sequence of viability kernels until time k satisfies 



V = V° and V fc+ i = {x £Y k j 3u € U , (x, u) G B and f(x, u) € Y k } 



(24) 



By (|23p . such an algorithm provides approximation from above of the viability kernel as follows: 



V(/,B) C Pi V fc = lim |V fc . 

1 1 k— >+oo 



(25) 



fceN 



Conditions ensuring that equality holds may be found in (jSaint-Pierrel . 11994! ) . Notice that, when the 
decreasing sequence (Yk)k&n of viability kernels up to time k is stationary, its limit is the viability 
kernel. Indeed, if Yk = Yk+i for some k, then Yk is a viability domain by (124(1 . Now, by ( I19|) . 
V(/, D) is the largest of viability domains. As a consequence, Yk = Y(f, B) since V(/, B) C Yk 
by (f23l) . We shall use this property in the following Sect. [Bj 



B Viable control of generic nonlinear ecosystem models with har- 
vesting 

For a generic ecosystem model ([1]), we provide an explicit description of the viability kernel. Then, 
we shall specify the results for predator-prey systems, in particular for discrete-time Lotka-Volterra 
models. 

The acceptable set B in (fT7|) is defined by minimal biomass levels By > 0, B\ > and minimal 
catch levels C\ > 0, C\ > 0: 

B = {(y,z,v,w) GR 4 \y > B y , z>B\, vy > Cj, wz > & z } . (26) 
B.l Expression of the viability kernel 

The following Proposition [5] gives an explicit description of the viability kernel, under some condi- 
tions on the minimal levels. 

Proposition 5 Assume that the function R y : M 3 — > K is continuously decreasing in the control 
v and satisfies lim„_s. +00 R y (y, z, v) < 0, and that R z : R 3 — > M is continuously decreasing in the 
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control variable w, and satisfies \\m w ^. +00 R z {y,z,w) < 0. If the minimal levels in (|26p are such 
that the following growth factors are greater than one 

R y {B y ,B\,% > 1 and R^Bl^i) > 1 , (27) 

the viability kernel associated with the dynamics f in ^ and the acceptable set D in (f26|) is given 
by 

Y(f,B) = l(y,z)\y>B y , z>B\, yR y (y, z,% > B\, zR z (y, z, > B\ 1 . (28) 



Proof. According to induction (|24l) . wc have: 



Vo 



{(y,z) y>B b z>B\ }, 



Vi 



(y, z ) 



V > By, z > B\ and, for some (v, w) > 0, 
vy>Clwz>ClyR y (y,z,v)>BlzR z {y,z,w)>Bl 



CI 



CI 



y>B> y ,z>B 9 z ,yR y {y,z,^)>Bl,zR z (y,z,^)>Bl 
" y " z 



because v M- R y {y, z,v) and w^R z {y,z,w) are decreasing 

and thus we may select v = — -, w = — . 

V z 

Denoting y' = yR y (y, z,v), z' = zR z (y, z,w), we obtain, 



(y, z ) 



y > B y , z>B\ and, for some (v, w) > 0, 



vy > C y , wz > CI 

y' > B y , y'R y {y', z ',^L) > B\, z> > B b z , z'R z (y',z',£) > B\ 
We shall now make use of the property, recalled in Sect. [3] that when the decreasing sequence (Vfe)fc 6 N of 
viability kernels up to time k is stationary, its limit is the viability kernel V(/, ID). Hence, it suffices to show 
that Vi C V 2 to obtain that V(/, D) = Vi. Let (y, z) € Vi, so that 

y>B b y , z>B\ and yR y (y, z, % > B\, zR z (y, z, Q) > B\ . 

y y y z 

Since R v : R 3 R is continuously decreasing in the control variable, with \im v ^. +00 R y (y, z,v) < 0, and 
since yR y (y,z, > B y , there exists a v > ^ (depending on y and z) such that y' = yR y (y, z,v) = B y . 
The same holds for R z : R 3 -> R and z' = zR z (y, z, w) = B\. By ([27). we deduce that 



y'Ryiy', z', ^) = B y Ry(By, B b z , ^) > B\ and z'R z (y', z< , Q) = B b z R z {B b y ,B\, J) > B\ 



The inclusion Vi C V2 follows. 
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Corollary 6 Suppose that the assumptions of Proposition are satisfied. Denoting 



v(y,z) = max{« > -f | yR y {y,z,v) = y b } 



v 



w(y,z) = max{iu > j zR z (y,z,w) = z b } , 
the set of viable controls is given by 



U V (/,D)(y>^) = S ( v , w ) 



v(y, z)>v>^, w(y, z)>w>^f, 



y>Ry(y>, z', ^)>y\ z'R z (y', z', ^) > z 



where y' = yR y (y,z,v), z' = zR z (y,z,w). 



B.2 Proof of Proposition [4] 

Proof. By (j9]), and the property that both R y and R z are decreasing in the control variable, the quantities 
exist. 

Also since both R y and R z are decreasing in the control variable, we obtain that 



Cl:* 



^77- >-Ry(B y ,Bl^L r ) = l and R z (B b y ,B b z , & **^ Zq) )> W b y , B b z ,^) 



Bl 



Cl 
B b 



y v 

To end up, the above inequalities and the assumption that yo > B y an( l z > £^ allow us to conclude, thanks 
to Proposition [51 that (yo,zo) belongs to the viability kernel V(/, D) given in (|28p . 

In other words, starting from the initial point (y(to), z(to)) = (yo,Zo), there exists an appropriate har- 
vesting path which can provide, for all time, at least the catches ©. □ 
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